Dislocations in the ground state of the solid-on-solid model on a disordered substrate 
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We investigate the effects of topological defects (dislocations) to the ground state of the solid-on- 
solid (SOS) model on a simple cubic disordered substrate utilizing the min-cost-flow algorithm from 
combinatorial optimization. The dislocations are found to destabilize and destroy the elastic phase, 
particularly when the defects are placed only in partially optimized positions. For multi defect 
pairs their density decreases exponentially with the vortex core energy. Their mean distance has a 
maximum depending on the vortex core energy and system size, which gives a fractal dimension of 
1.27 ± 0.02. The maximal mean distances correspond to special vortex core energies for which the 
scaling behavior of the density of dislocations change from a pure exponential decay to a stretched 
one. Furthermore, an extra introduced vortex pair is screened due to the disorder-induced defects 
and its energy is linear in the vortex core energy. 



I. INTRODUCTION 

At low temperatures the physics of crystal surfaces on 
disordered substrates is dominated by the randomness 
rather than thermal fluctuations. In 2+1-dimcnsions 
this elastic surface is expected to have a roughening 
phase transition at a critical temperature T c from a ther- 
mally roughs-phase for T > T c to a superrough phase 
for T < T c ETE1, corresponding to a height-height cor- 
relation function log(r) and log 2 (r) respectively. The 
log 2 (r)-superrough behavior was numerically confirmed 
at finite temperature via Monte Carlo simulationso as 
well as in the limit of a vanishing temperature via exact 
ground state|-calculation using combinatorial optimiza- 
tion methodaJLl. 

In this paper, we study the stability of the low- 
temperature (glassy) phase of thej-|Solid-on-solid model 
(SOS) on a disordered substratetXEl with respect to the 
formation of topological point-like defects. We also con- 
sider the density of defects and the screening effect of 
pre-existing pairs to an introduced extra pair and allow 
for a vortex-core energy. 

The SOS model on a disordered substrate is given by 
the uniformly distributed substrate height di £ [0, 1] and 
the integer crystal height rij on a simple cubic L x L- 
lattice G with periodic b.c. and lattice site i as schemat- 
ically shown in Fig.|l|. The hi — n, + di denotes the total 
surface height at site i and the SOS model Hamiltonian 
is defined by 



(ij) 



(1) 



where the sum runs over all nearest-neighbor pairs (ij). 
To calculate the ground state of the SOS Hamiltonian (fil) 
we introduce the crystal height-differences n*j = rii — n,j 
(integer) and substrate height-differences d*- = dj — di 
(e [—!,+!]) along the links k = (i,j) on the dual lattice 



^ -x,y 

FIG. 1. Height profile hi = m + di in the (2+1) SOS 
model, where di £ [0, 1] are the random offsets of the disor- 
dered substrate and Hi the crystal heights (arbitrary integers) 
on the lattice site i. 



G* . Thus we get the following SOS Hamiltonian for the 
dual space 

ff({n£}) = 5>£-d£) a . (2) 



The minimal (optimal) energy configuration {n^} m i n 
will just be the closest integer n£ to d£ for all links 
& = (i-j)- On the other hand, since the n£ describe 
height-differences in the scalar field given by the rij their 
sum along any oriented cycle on the surface around site 
i has to be zero, i.e. the lattice divergence of n* has to 
vanish for each site i: 



(V-n* 



0. 



(3) 



Note that rii can be considered to be a potential and 
n|,- as its force field. Obviously, for a typical disordered 
substrate the minimal configuration {n£} m i„ violates the 
mass balance constraint (^). Fig.^ shows an example of 
a disordered substrate with substrate height di = 0.0, 
0.2, 0.4 and 0.6. Consider the differences d£: across the 



dashed line we have d£ = 0.6 and 



d£| < 0.5 elsewhere. 



Consequently, the absolute minimum-energy configura- 
tion without any balance constraint is given by n£ = 1 
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FIG. 2. Example of a disordered substrate heights di in a 
random-surface model with a single dislocation pair connected 
along a straight line of size L (dashed line). The optimal sur- 
face without dislocations would be flat, i.e. m = for all sites 
i, however, allowing dislocations would decrease the ground 
state energy (see text). 

and n£ = respectively. With respect to the balance 
constraint ([?]) the only feasible optimal solution (ground 
state) is a flat surface, i.e. n£ = for all links k = (w). 
On the other hand, dislocations of Burgers chargeEj b 
can be introduced if one treats the height field hi as a 
multi-valued function which may jump by b along liues 
that connect two point defects (i.e. a dislocation pair)lill. 
Therefore, for the given example (FigJ|) it should be clear 
that the minimal configuration {n* k } rn i n (see above) is ex- 
actly the optimal (i.e. ground state) configuration with 
one dislocation pair. One of the two defects has a Burgers 
charge b = +1 and the other one b = — 1. The pair is con- 
nected by a dislocation line (dashed line in Fig. ||) along 
which one has n* = 1. This already demonstrates that 
due to the disorder the presence of dislocations decreases 
the ground state energy and .a proliferation of defects ap- 
pears. Alternatively, in Ref.E a dislocation pair (excited 
step) was introduced by fixing the boundary to zero and 
one. 



II. DEFECT PAIRS IN THE SOS MODEL 

The defect pairs in the disordered SOS model are 
source and sink nodes of strength Arb and —b, respec- 
tively, for the network flow field riicHa, which otherwise 
fulfills (V • n*)j = 0, i.e. we have to modify the mass 
balance constraint (||) as follows 

, *\ _ ) , no dislocation at i 

^ ' % ~ 1 ±6 , dislocation at i ^ ' 

Thus the ground state problem is to minimize the Hamil- 
tonian (|J) subjected to the mass balance constraint 
(Q) which-can be solved by the successive-shortest-path 
algorithnflla'Ll. In the following we concentrate on defect 
pairs with b = ±1. 

The defect energy AE is the difference of the mini- 
mal energy configuration with and without dislocations 
for each disorder realization, i.e. AE = E\ — Eq. More 




FIG. 3. Graph of a L x //-lattice with periodic b.c. for 
the implementation (a) of one fixed defect pair and (b) of a 
partially optimized pair. Both are separated by L/2. The en- 
ergetic costs are Cfc(n^) = (n£— d£) 2 at the dual site k = (i,j). 
Dislocations are induced by two extra nodes s and t, which 
are connected with the possible positions of the defects (big 
dots). 

precisely, for the configuration with N defect pairs of 
Burgers charge b = ±1 we introduce two extra nodes 
s and t with n s = +N and n t = —N respectively and 
connect them via external edges or bonds with particular 
sites of the lattice depending on the degree of optimiza- 
tion: (a) with two sites separated by L/2 (Fig.||(a)), (b) 
the source node with one site i and the sink node with the 
sites on a circle of radius L/2 around i (Fig.||(b)) and (c) 
both nodes with the whole lattice. Case (a) corresponds 
to a fixed defect pair, (b) to a partially optimized pair 
along a circle, both separated by a distance L/2, and (c) 
to a completely optimized pair with an arbitrary separa- 
tion. In all cases the energy costs for flow along these 
external edges are set to a positive value in order to en- 
sure the algorithm to find the optimal defect pair on the 
chosen sites. These "costs" have no contribution to the 
ground state energy. In case of multi pairs we always 
use graph (c). Here, the optimal number N of defects 
in the system is gradually determined starting with one 
pair (N = 1) with a vortex core energy 2E C and check- 
ing whether there is an energy gain or not. If yes, add 
a further pair (with 2E C ) and repeat the procedure until 
there is no energy gain from the difference of the ground 
state energy between two iterations. 



III. SINGLE DEFECT PAIR (JV = 1) 



We study the defect energy AE and its probability dis- 
tribution P{AE) on&LxL lattice with L = 6, 12, 24, 
48, 96 and 192 and 2 ■ 10 3 — 10 5 samples for each size and 
consider the three cases (a)-(c) (see above). With an in- 
creasing degree of optimization a negative defect energy 
AE becomes more probable and its probability distribu- 
tion P(AE) differs more and more from the Gaussian 
fit, Fig.[|. The resulting disorder averaged defect energy 
[AE]di S scales like 
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FIG. 4. Probability distribution P(AE) of a large-scale topological excitation with a Gaussian fit for different optimizations: 
(a) for a fixed defect pair, (b) for a partially optimized pair and (c) for a completely optimized pair with different system sizes L. 



ln(L) , fixed defect pair 

[AE] dls ~ <( -0.27(7) • ln 3/2 (L), partially optimized 

-0.73(8) • ln 3/2 (L), completely optimized 



and its related variance a like 

ln(L) , fixed defect pair 
a(AE) ~ ^ ln 2/3 (L), partially optimized 

ln 1 ' 2 ^), completely optimized 



(5) 



(6) 



where the exponents are approximations for the best data 
collapse. The defect energy indicates that for the opti- 
mized cases dislocations can proliferate due to thermal 
fluctuations and melt the elastic superrough phase. Fur- 
thermore, for a growing degree of optimization the scaling 
amplitude of [Ai?]^ increases. 

The mean length (mass) Idl of the line connecting the 
two defects scales with the system size L according to the 
fractal dimension 



d f = 1.28 ±0.02 

for the fixed and partially optimized situation. 



(7) 




FIG. 5. Finite-size-scaling relation of the probability dis- 
tribution P(Idl) of the mean distance Idl between two op- 
timally placed dislocations for system size L — 6, 12, 24, 48, 
96 and 192. The data collapse for P L (l DL ) ~ 1/L ■ p{l DL /L). 



For the completely optimized case Fig j5] shows a proba- 
bility distribution P(Idl), which behaves like 



Pl{1 



DL 



1 

L 



Idl 
17 



(8) 



IV. MULTI DEFECT PAIRS (N > 1) 



Next, we study the effect of a uniformly given vortex- 
core energy E c to the system of multi deject pairs (N > 1) 
as a simplification of the real situation with a distribu- 
tion of E c . As shown in Fig.|fj|(a), the density p of defects 
decays exponentially with an increasing E c , i.e. 



P(E C 



^-E c /E )° 



(9) 



For the Eq and a we can distinguish between two inter- 
vals of E c which refers to a stretched and a pure expo- 
nential decay, respectively. In detail we have 



E c e 


E 


a 


[0,oo [ 


0.75 ±0.15 


0.6 ±0.2 


[0,E™ ax (L)[ 


0.45 ±0.03 


1 



The upper limit E c max (L) corresponds to the maximal 
mean length Idl for each system size L, c.f. Fig.[|(a) and 
(b), and scales like 



E r 



(const. + 0.47 ±0.02) ■ ln(L)) 



3/2 



(10) 



Moreover, we found the same scaling behavior for the 
vanishing defect energy, i.e. [A-E]^ s = 0: 



E c0 ks (const. + (0.47 ± 0.01) • ln(L)) 



3/2 



(11) 



From the plot of the maximal mean length Idl (Fig.^(b)) 
vs. the system size L, i.e. Idl{E^ ii1x ) ~ L d f , the fractal 
dimension df is given by 



d f = 1.267 ±0.07, 



(12) 
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FIG. 6. (a) Density p of defects with respect to the vortex 
core energy E c for different system sizes L = 6 - 48 and 10 3 
up to 10 4 samples. The log-lin plot indicates an exponential 
decay of p. Simultaneously, (b) the mean distance Idl of all 
dislocation pairs vs. the vortex core energy E c . Comparing 
(a) and (b) one sees that the maximal lengths Idl occurs at 
the cross-over energy E™ ax (see text). 

close to the one of the above single line situation (df s» 
1.28). 

Finally, we focus on the effect of introducing an extra 
fixed defect pair separated by L/2 to an already (com- 
pletely) optimized configuration with a vortex core en- 
ergy E c . This extra pair costs 

AE fix =E[ + 2E c -E 1 , (13) 

where E\ denotes the ground state energy for N (pre- 
existing optimal) pairs and E[ the energy for N + 1 op- 
timally placed pairs, both for the same disorder configu- 
ration {di}. As plotted in Fig.|?], AEfi x is constant in L, 
but linear in E c , i.e. 

AE fix (L) = (0.17 ± 0.02) + (4.35 ± 0.02) • E c . (14) 



FIG. 7. Defect energy AE of a single (introduced) defect 
pair vs. the system size L (L = 6, 12, 24, 48, 96) in a system 
with and without an already optimal number of dislocations 
for different vortex core energies E c . 

Thus, one obtains a screening effect of the defect- 
defect interaction due to disorder- induced dislocations. 
In comparison, Fig.|^ also shows the case for a single pair 
(N = 1) without pre-existing pairs as studied in section 
III. 



V. RELATED MODELS 

A similar picture of the effect of dislocations to a ran- 
domly pinned elastic media at T = were found for_Dther 
discrete models, theJully-packed loop (FPL) mode£3 and 
the matching modelE-3, both on a bipartite hexagonal lat- 
tices with a linear energetic cost function and periodic 
b.c. 

In the case of a single fixed defect pair we found the 
same ln(L)-behavior x>£-jthe defect energy as for the ex- 
citation step in RefE'EJ, but got a smaller fractal di- 
mension d f = 1.28(2) rather than d f = 1.35(2)0. The 
disorder-induced dislocations turned out to destroy the 
quasi-long-range order of the elastic phase due to a neg- 
ative scaling behavior of defect energy AE with respect 
to L for optimally placed defects, i.e. AE ~ — ln 3 ^ 2 (Ll, 
in good agreement with the results of the FPL modelt3. 
When taking into account screening and an uniform vor- 
tex core energy E c in addition to the energy balance AE, 
one finds that the energetic costs AEfi X of an introduced 
fixed pair does,not depend on the system size L, as shown 
for the FPLE2I and matchingta model. In addition we 
found for the disordered SOS model a linear dependence 
of AEfi x on the vortex core energy E c (Fig. [?]). One 
concludes that this extra pair is screened by the pre- 
existing defect pairs. For the exponentially decay of den- 
sity p of the dislocation pairs we distinguished between 
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(a) the whole range of the vortex core energy E c and (b) 
a range with an upper limit E™ ax {L), for which (latter 
case) the mean defect length Idl was found to be maxi- 
mal. The case (a) corresponds to a = 1 and Eq = 0.45(3) 
and the maximal length / dl related to the cross-over en- 
ergy E™ ax (L) behaves as Idl ~ L f with the fractal 
dimension df — 1.267(7). Both results were also found 
in Refill. For the case (b) we get a w 0.75iclose to 2/3) 
and Eq » 0.6 in good agreement with Ref.t! 

Finally, we relate the SOS model (0) to the continuum 
description of a randomly pinned clastic medium on large 
length scales given by the sine- Gordon model Hamilto- 



H= / d 2 r 



f(V«(r)) 2 



w cos (27r(«(r) — d(r))) 



(15) 



where K is the elastic constant and d(r) a random field 
out of [0,1]. The first term represents the elastic en- 
ergy E e i and the second one the random pinning energy 
Epi n . The model is known to describe a weakly disturbed 
vortex lattice in a thin two-dimensional (2IX) pSupercon- 
ducting film introduced by a parallel fieldll3~bZl. XJther 
experimental realizations are charge, density wavesES and 
Wigner crystallization of electrons!!! 

The relation to the SOS model is as follows: in the 
limit of an infinite coupling strength w — > oo and T = 
the sine-Gordon model maps onto a lattice SOS model 
Eq.fll]), as the cosine-term of Eq.(p"5|) fojapsn-the displace- 
ment field u(r) to be u(r) — d(r) + n(r)pH'EI ) where n(r) 
is an integer. One can identify the u(r) as the continuous 
height field h(r) of the SOS model. 

Th e r . esYilt j3 of the analytical study of the sine- Gordon 
modelE3'EHIc3 are in good agreement with our results, but 
only refer to the cases of fixed and completely optimized 
pairs. Furthermore, these studies allow another inter- 
pretation of the defect energy AE and .density p. From 
the calculatiop_pf the elastic energy EjEi and defect en- 
ergy A£cj'E3'E3 one gets that for a fixed pair the elas- 
tic energy E e i dominates the pinning energy E p i n , i.e. 
AE ~ E e i, and for the completely optimized pair the 
situation is vice versa, i.e. AE ~ E P i n . The result- 
ing scaling behavior is found to be AE ~ ln(L) and 
AE ~ — ln 3 / 2 (L), respectivelyE3. The scaling behav- 
ior of the fixed dislocation pair in presence of pinning 
disorder is essentially equivalent to the one of a fixed de- 
fect pair at finite temperatures T without disorder, i.e. 
A£~ln(L) ~E$™(T). 

The density p of defects can be related to the 
length scale £d beyond which the dislocations become 
unpaired 2 !] since for a « 0.6 the density p, Eq. (||), shows 
the same scaling behavior as £d in the case of low tem- 
peratures and large core energy E c , i.e. E c 3> K In(fD). 
For E c ks we found large densities p and one is prob- 
ably out of the regime given by E c ^> KIii^d)- This 
would possibly explain the occurrence of the stretched 
exponential behavior close to E c « as seen in Fig.||(a). 



To summarize we have studied the effect of disloca- 
tion pairs on the ground state properties of the SOS 
model on a disordered substrate. For a fixed position 
of the dislocation pair a distance L apart we found that 
on average t he defect costs an energy proportional to 
InL, in agreement with the findings for the energy costs 
for a steplike_excitation step with fixed end points re- 
ported earlieru and also in agreement with-jjecent results 
for other two-dimensional lattice modelaia'cj. On the 
other hand if we optimize the position of the dislocation 
pair we showed that it gains energy, namely an amount 
proportional to In L with an exponent tp around 3/2 as 
predicted tw-,scaling arguments and also observed in the 
FPL model! When introducing a penalty for the topo- 
logical defects (i.e. a core energy) we showed that the 
density of defects vanishes exponentially as a function of 
this core energy, which is in agreement with the results 
for the FPL model! Finally we also demonstrated that 
a dislocation pair is screened by the presence of other 
dislocations in the system. 
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